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ANALYSIS OF A REDUCED-ORDER HDG METHOD EOR THE 

STOKES EQUATIONS 


ISSEI OIKAWA 


Abstract. In this paper, we analyze a hybridized discontinuous Galerkin(HDG) 
method with reduced stabilization for the Stokes equations. The reduced sta¬ 
bilization enables us to reduce the number of facet unknowns and improve the 
computational efficiency of the method. We provide optimal error estimates in 
an energy and norms. It is shown that the reduced method with the lowest- 
order approximation is closely related to the nonconforming Crouzeix-Raviart 
finite element method. We also prove that the solution of the reduced method 
converges to the nonconforming Gauss-Legendre finite element solution as a 
stabilization parameter r tends to infinity and that the convergence rate is 
O(r-l). 

discontinuous Galerkin method and hybridization and Gauss-Legendre element 
and Stokes equations 


1. Introduction 

The aim of this paper is to propose and analyze a reduced-order hybridized dis¬ 
continuous Galerkin(HDG) method for the Stokes equations with no-slip boundary 
condition: 

—Am -|- Vp — f in fl, 

(1.1) divM = 0 in 11, 

M = 0 on d^l, 

where ft C W^{d = 2,3) is a convex polygonal or polyhedral domain and / € 
:= is a given function. For the Stokes problem, various HDG 

methods were already proposed and studied [HlinillslEllinillolEllIllEIlIIH]. 
We also refer to m for an overview. The method we investigate in this paper is 
the HDG-IP method proposed by Egger and Waluga in [21] . The HDG-IP method 
is based on the gradient-velocity-pressure formulation of the Stokes equations. We 
remark that the HDG method of the local discontinuous Galerkin type (22] [12] is 
close to the HDG-IP with a slight difference of a numerical flux. 

A reduced stabilization was introduced to the DG methods and was analyzed for 
elliptic problems [6] |7]. In [4], Becker et al. studied the reduced-order DG method 
for the Stokes equations and analyzed the limit case as a stabilization parameter 
tends to infinity. In HD, Lehrenfeld proposed a reduced-order HDG method for the 
Poisson equation and the Stokes equations. Lehrenfeld also remarked that the con¬ 
vergence rate of the method is optimal, however, error analysis was not presented. 
In [23], for the Poisson problem, the author provided the optimal error estimates 
and showed the reduced-order HDG method with the lowest-order approximation is 
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closely related to the nonconforming Crouzeix-Raviart finite element method. Re¬ 
cently, Qiu and Shi analyzed the reduced methods for linear elasticity problems [25] 
and convection-diffusion equations [23]. 

The reduced-order HDG method uses polynomials of degree fc -|- 1 and k to ap¬ 
proximate element and hybrid unknowns, respectively, whereas the standard HDG 
method uses polynomials of degree k for both unknowns. Although both the meth¬ 
ods have the same number of globally coupled degrees of freedom, the reduced 
method can provide higher order convergence than the standard method. This is 
the main advantage of the reduced-order HDG method. For the Stokes problem, 
when we use polygonal or polyhedral elements, the reduced method is indeed better 
than the standard method in term of convergence orders. The standard method 
using polynomials of degree k for all unknowns obtains the suboptimal convergence; 
the orders are A:-I-1 for velocity without postprocessing and k+1/2 for the gradient 
and pressure, according to [13]. In contrast, the reduced method uses polynomials 
of degree fc -|- 1 for velocity and polynomials of degree fc for the hybrid part of 
velocity and pressure and can achieve the optimal order convergence. 

In this paper, we provide optimal error estimates of the reduced method for 
the Stokes problem. Since we need to use a weaker energy norm in our analysis, 
it is necessary to modify the error analysis of the standard method. We note 
that the main difficulties can be overcome by the techniques used in the author’s 
previous work [23] and the discrete inf-sup condition proved by Egger and Waluga 
[18] . We also show a relation between the reduced method and the Gauss-Legendre 
element (see [5] [26] for example). It is proved that the hybrid part of velocity and the 
pressure of the reduced-order HDG method with the lowest-order approximation 
coincides with those of the nonconforming Grouzeix-Raviart finite element solution. 
In the limit case as the stabilization parameter r tends to infinity, the solution of the 
reduced method converges to that of the nonconforming Gauss-Legendre method. 
The convergence rate is estimated to be This result is inspired by [d] 

Theorem 3], however, our proof is completely different and novel. 

The rest of this paper is organized as follows. Section 2 is devoted to the prelim¬ 
inaries. In Section 3, we introduce a reduced stabilization and present a reduced 
HDG method. In Section 4, we provide a priori error estimates in an energy and 
norms. In Section 5, some relations between the nonconforming Gauss-Legendre 
finite element method and the reduced method are shown. In Section 6, numerical 
results are presented to confirm our theoretical results. 


2. Preliminaries 

2.1. Meshes and function spaces. Let {Th}h be a family of shape-regular trian¬ 
gulations of D and define F/i = UiceTh edges in Th- The 

mesh size of Tk is denoted by fc, namely fc := maxKeTh ^k, where hx = diamAT. 
The length of an edge e € Sh is denoted by he- 

We use the usual Lebesgue and Sobolev spaces; L^(D), L‘^{Th) and ifc™(Sl), and 
also Lo(^) = {? S : J^gcix = 0}. We introduce piecewise Sobolev spaces 

H^{Th) = {u G L^(D) : v\k G H'^{K) VAT G Th}- For vector-valued function 
spaces, we write them in bold, such as L'^{Vl) = and 

The usual inner product is denoted by (•, ■)q- Let us define the piecewise inner 
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products by 

{u,v)rn= X! / X! / 

K^Th KGTh 

Let PkiTh) and Pk{£h) denote the space of element-wise and edge-wise polynomials 
of degree k, respectively. We employ = Pk+i{Th), — Pk{£h) H {v S 

L'^iTh) : v = 0 on 517} and = Pk{Th) H Lq{£1) as finite element spaces, which 
we call Pk+i-Pk/Pk approximation. The L^-projection from niCGT;. L‘^{dK) onto 
Wx^rh Pk{dK) is denoted by P^, and I stands for the identity operator. 

2.2. Norms and seminorms. As usual, we use the Sobolev norms |i?|m = |'i’|iT™(n) 
and IliilIm.D = ||i’||ffm(n) for a domain D. The L^-norm is denoted by ||i;|| = 
||'w||o,D = ||i’||l 2 (£)). The energy norms are defined as follows: for {v,v) G H^{Th)x 
LH^h), 

\iiv,v)r 

\liv.v)\ll 
\Uv,v)\ll 

where 

KGTh eCdK 

\v\Im = X 

KgU 

\'^\2M = X 

Ken 

The symbol r is a stabilization parameter which will be defined in Section [3l The 
parameter-free energy norms ||| • ||| and ||| • |||/i are used to analyze the convergence rate 
with respect to the mesh size h. We need the parameter-dependent energy norms 
in order to analyze the proposed method when r ^ oo. We also use the stronger 
norm 

Ikll^ = IkiP + X ^Kkli.R-- 

Ken 

By the inverse inequality [I], we see that the two energy norms are equivalent to 
each other on x Vjf, i.e., 

(2.1) |||(ii.,^l.)||U< 11(11.,ll.)|||<C|||K,i;.)||U, 

(2.2) \liVh,Vh)\lh,T < l|(ll.,v.)|lr < C\l{Vh,Vh)\lh.r 

for some constant C > 0 independent of h. Similarly, it holds that < ||g.||. < 
Cllg^ll for all Qh G Qk- In the following, the symbol C will stand for a generic 
constant independent of the mesh size h and the stabilization parameter r. 


= bll,. + 1^12,. + l(^’,«)lf, 

= \v\lh + \v\ih + l(^,?^)lj- r, 

= bll,. + l(^,«)lj^ 

= \v\lh + l(«,«)lj- T. 
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2.3. Approximation property. The approximation property in the energy norm 
holds as well as in the standard energy norm. 

Theorem 1. If tp € and tt G then we have 

(2.3) inf ^ \l{ip - Vh,ip\rh-Vh)\l < 

(2.4) inf IItt - g^ll < C'/i'=+VlH'=+i(n)- 

qhdQl 

Proof. We refer to [53]. □ 

3. A REDUCED-ORDER HDG METHOD 

In this section, we present a reduced-order HDG method based on the HDG 
method proposed by Egger and Waluga in m- By taking the L^-projection onto 
the polynomial space of lower degree by one in the stabilization term of the standard 
method, we obtain the reduced-order HDG method: find {uh,Uh,Ph) G x 

X Q'l such that 

(3.1a) ah{uh,Uh;Vh,Vh) + bh{vh,Vh;Ph) = {f,Vh)n G x 

(3.1b) bh{uh,Uh;qh) =0 'iqh & Qt, 

where the bilinear forms are given by 

'^hi '^h') ^'^h)Th '^h '^h)dTh '^h '^h}dTh 

+ {Thf^Pk{Uh - Uh), Pk{Vh - Vh))dTh, 
bhivh,Vh-,Ph) = -{divVh,Ph)Th - i^h - Vh,Phn)dTh- 

Here r is a stabilization parameter assnmed to be greater than or equal to one 
and sufficiently large. We recall that P^, which is defined in Section 2.1., is the 
L^-projection onto the edge-wise polynomial space of degree k. 

Remark 2. In the two-dimensional case, we can easily implement the redueed 
method by using a reduced-order quadrature formula in the eomputations of the 
reduced stabilization term, see [531 Lemma 5] for details. 

4. Error analysis 

In this section, we provide the optimal error estimates of the method in both 
the energy and norms. To do that, we first show the consistency of the method, 
the boundedness of Oh and bh, and the coercivity of o^,. In addition, the discrete 
inf-sup condition of b^ is proved based on the results of |T3] . 

4.1. Consistency. We state the consistency and adjoint consistency of the method. 

Theorem 3. Let {u,p) be the exact solution of the Stokes equations (11.11) . Then 
we have 

(4 1) ah{u,u\Y,^-,Vh,Vh) + hh{vh,Vh]p) = {f,Vh)Q. V(uft,,u,,) G x , 

bh{u,u\r^-,qh) =0 'iqh&Qh- 


Proof. Since it — M|r;, = 0 on T^, we can easily see that the consistency (14.11) 
holds. □ 
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Let {uh,Uh,Ph) G X X be the solution of the method (13.ip . From 

the consistency, the Galerkin orthogonality follows immediately: 

(4.2) 

ah{u - Uh,u\r^ - Uh]Vh,Vh) + bh{vh,Vh;p - ph) = 0 \/{vh,Vh) G x 

'^\rh Qh) — 0 G Q 

Due to the symmetricity of ah-, we readily see that the adjoint consistency also 
holds: 

ah{vh,Vh;u,u\r,,)+bh{vh,Vh;p) = if,Vh)Q 'i{vh,Vh) & x , 

bh{u,u\r^-,qh) =0 \/qh & Qh- 

4.2. Boundedness and coercivity. We first prove the boundedness of an and 

bh- 

Theorem 4. Let (|, ^) = {w + Wh, + Wh) and {rj, fj) = {v + Vh, u|r;. + Vh), 

where {wh,Wh), {vh,Vh) G x and w,v G Hq{LI). Then there exists a 

constant C > 0 independent of h and t such that 

(4.4) |a.(|,|';r,,r^)|<Cr|||(4,|')||||||(r,,i7)|||. 

With respect to the parameter-dependent energy norm, we have 

(4.5) |a.(l,f;r/,77)|<C|||(^,allWII(r7,^7)|||.. 

Proof. By the Schwarz inequality, the first term of ah is bounded as 

(4.6) |(V|,Vr7)rJ < llliAlr/li./,. 

We estimate the second term. Note that 

{9r,^,v-v)dn = {dn^^Pkifj - v))dn + (^«4,(i - Pk){f]-r]))on 
= {9-n^, Pkiv - 'n))dTh - {9n^, (I - Pk)Vh)aTh- 
Since {d-nW, (I - Pk)v)9rh = {d-nWh, (I - Pk)v)dTh = 0: if follows that 

{d-n^, (I - Pk)v)drh = 0. 

Using this, we deduce that (I - Pk)vh)dTh = {dn$, (I - Pk)ri)dTh- Then we 

have 

(4.8) 

\{dn^,V -v)dn \ = \{dn$.,Pk{Vh-Vh) - (I - Pk)v)dn\ 

< CmaxlGr-i/^KI^I? ;, + /i2|^|2_^)i/2(|(r7,77)|2^ -h 

<cm,mr\Kv.mr, 

where we have used the trace inequality and the following estimate (see [53] for the 
proof) 

(4.9) \iv,v)\i <C\v\ih- 

The stabilization term is bounded as 

(4.10) \{Thf^Pk{i- ^),Pkifj - 'n))arj < 

From (|4.7I) . (14.81) and (14.101) . we obtain the boundedness (14.5F With a slight mod¬ 
ification, we can also show that (iTl) . □ 
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Theorem 5 (Boundedness of bh). Let (^,0 be the same as in Lemma and 
r = q + qh with q G H^{Th) H Lq{LI) and qh G Then there exists a constant 
C > 0 independent h and r such that 

(4.11) |fe.(tl';r)|<C|||(e|')|||||r|U. 

In particular, in the case of r = qh G Q^, we have 

(4.12) l^>/.(€,l';9.)l<C'|||(C,l')llllk.||. 

Proof. By the Schwarz inequality, we have 

(4.13) l(div^,0rj < |l|i,/«lkl|. 

Next, we estimate the second term of the bilinear form bh- In a similar manner of 
(EH), we have 

{^-i,rn)orh = (Pfc(f- i),rn)Qrh + ((I ^ Pfe)(^ “ Oyrn)gr„ 

(4.14) 

= (Pfc(€ - i),rn)drH - ((I - Pk)wh,rn)drH- 

Since ((I - Pk)w,rn)Qrh = 0; we get 

((I - Pk)wh,rn)srh = ((I - Pk){w PWh),rn)ar,,, = ((I - Pk)$.,rn)Qrh- 

H6I1C6 

\(i-i:rn)grh \ = KPfcd'- I) - (I - Pk)^:rn)grJ 
<cm,i)\! + K\lhy^^r\\h, 

where we have used the trace inequality and (14.91) . Consequently, we obtain the 
inequality (I4.11I1 . From the inverse inequality, (I4.12|) follows immediately. 

□ 

Theorem 6 (Coercivity). Assume that r is sufficiently large. There exists a con¬ 
stant C > 0 independent of h and r such that 

(4.15) ah{Vh,Vh;Vh,Vh) > C'|||(U;,,U;,)|||^ V{Vh,Vh) e X Vh- 

In particular, we have 

(4.16) ah{vh,vh;vh,vh) > C\l{vh,vh)\f V{vh,vh) G x . 

Proof. We note that 

(4.17) {dnVh,Vh - Vh)aTh = {9nVh, Pkivh - Vhfau- 
Then it follows that 

(4.18) ah{Vh,Vh\Vh,Vh) > \Vh\lM “ ‘^\{dnVh, Pk{Vh - Vh))dTh\ + \{'^h,Vh)\lr- 
By the trace and inverse inequalities and Young’s inequality, we have 

(4.19) 2\{dnVh,Pkivh - Vh))aTh \ < C {e\vh\l^h + s~^T-^\{vh,Vh)\lr) 
for any e > 0. From (14.1811 and (I4.19|l . it follows that 

(4.20) ah{vh,Vh]Vh,Vh) > (1 - Ce)\vh\‘i^h + (1 “ Ce~^T~^)\{vh,Vh)\lr- 

We can take e = and deduce that, by assuming r > 4(7^, 

ahivh,vh;vh,vh) > (l - \l{vh,Vh)\ll,r > ^\livh,Vh)\ll^r- 

By the inverse inequality, we have |||(u/j,u/j)|||,- < C'\l{vh,Vh)\lh,T for some positive 
constant C", which completes the proof. □ 
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4.3. The discrete inf-sup condition. To prove the discrete inf-sup condition of 
the bilinear form bh, we introduce a Fortin operator. The main idea and techniques 
for constructing the Fortin operator are entirely based on [TB]. The global L^- 
projection operators 11^ : —>• Vj!^ and 11^ : —>• are defined by 

{n’^^v)\K = u>k{v\K) for K e %, 

(n^u)l, = n^(n|e) force 

where 11^ and Ilg are the L^-projections onto Pk{K) and Pk{e), respectively. We 
define the Fortin operator by 

In the following, we show this operator satisfies the Fortin properties. 


Theorem 7. For all v e iTg(n), we have 

(4.21) qn) = -(divu, qh)n ^qn e Q^. 

Moreover, there exists a constant C > 0 such that 

(4.22) |||(n^+ii;,n^«)||U < c\v\,.h \/v e H^m. 

Proof. First, we prove (14.2111 . Using the Green formula and the property of 
projection, we have 

bk{Ul+^v, qn) = -(divn^+'i;, g,)r, - {U^v - U^+^v, qkn)an 

(4 23) = Vg,)r. - g,n)ar. 

= {v,Vqh)TH - {v,qhn)dTh 
= -{(lwv,qh)n- 


Next, we prove (I4.22|) . Note that \n.^^v\i^h < and 


(4.24) 


h-l/2||p,(nfc„ _ n^+li;)||o.e < h-^/^\\v - Ul+^v\\o,e 

<C\v\i^h- 

From (14.231) and (j4.24|l . it follows that 

(4.25) I\inl+^v,n’f,v)i\l = \nl+^v\l^ + \{n’i+^v,nlv)\^ < c\v\l^, 
which completes the proof. 


□ 


By using the above results, we can prove the discrete inf-sup condition for the 
bilinear form bh- 


Theorem 8 (Discrete inf-sup condition). There exists a constant /3 > 0 indepen¬ 
dent of h such that 


(4.26) 


sup 


bh{vh,Vh;qh) 


> PhhW 


yqh e Ql 


Proof. It is well-known that the continuous inf-sup condition holds: there exists 
/?' > 0 such that 


sup 

^eJTj(a) 


(divt>,g)n 


> 


Vq e To(^)- 


(4.27) 
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For all g/t G Qft C Lq{Q), we have 


sup 

{vh,Vh)GV^+^ 


xV> 


bhivh,Vh;qh) ^ bhiUl+^v,U^v;qh) 


= sup 


•-h 

{divv,qh)n 

.eHiin) |||(n^+^,n»||U 


> C sup 

v£H^(Q) 

>l3'\\qHl 


{divv, qh)n 

\v\i,h 


The proof is complete. 


□ 


4.4. A priori error estimates. We prove optimal error estimates by using the 
results in the previous section. In this section, the stabilization parameter r is fixed 
to be a sufficiently large value. 


Theorem 9 (Energy-norm error estimate). Let {u,p) G HQ{Lt) x Lq{LI) be the 
exact solution of the Stokes equations (11.11) . and let (tt/j, Uh^Ph) G x x 
be the solution of the method dSH). Then we have 


(4.28) 


\l{u-uh,u\r^ - w?.)lll + \\P-Ph\\ 


< C 


( inf 



'^h,u\r^ 


^/i)lll+ infjb 



If{u,p) G ff'=+2(n) X £r'=+i(fl), we obtain 

(4.29) |||(M-M/„M|r^ -M/«)lll + \\P-Ph\\ < Ch'^+^ (|M|j/'=+2(n) + lbllff'=+i(n)) ■ 

Proof. Let Vh G , Vh G and rn G be arbitrary, and set r/h = Uh — Vh, 

f}h = Uh- Vh and 6h = Ph - qh- Then we have 

(4.30) 

ah{r]h,f]h;wh,wh) + bh{wh,Wh\5h) = F{wh,Wh) y{wh,Wh) G x V)f, 

bhiVh, fjh] Th) = G{rh) ^Th G Qh, 


where F : x M and G : Q'f —> M are defined by 

F{wh, Wh) = ah{u - Vh, ttlr^ - Vh;Wh, Wh) + bh{wh, Wh;P - qh), 

G{rh) = bh{u - Vh,u\ry^ - Uh, Vh)- 


By the boundedness of Oh and bh, we can estimate the dual norms of F and G as 
follows: 


||F|| = sup 

{wh,vih)£Vf^^ 


\F{Wh,Wh)\ 


\li'^h,Wh)\l 

< G{\l{u-Vh,u\r^ -«?.)lll + lb-'??.|l 


11*^11= sup < C\l{u - Vh,u\r^ - %)|||. 

rh&Qi 


(4.31) 
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By the coercivity (I4.15|l . we have 

c\l{vh,fih)\f < ahi'nh,fih;vh,vh) 
(4.32) = F{qh, f]h) - G{Sh) 

By the discrete inf-sup condition, we deduce that 

b[wh,Wh]5h) 


P\\Sh\\ < sup 

{■wh,wh)ev^~^^xv^ 




(4.33) 


< sup 

{■Wh.,iSh}(^V^^^xV^ 


F{wh, Wh) - ahi'nh,f]h; Wh, Wh) 

III Ill 


<||i^||+C'|||(77.,^7/.)|||. 

From (14.321) and (I4.33p . it follows that 

\l{Vh,fih)l\ + ||(5,i|| < C {l\{u - Vh,u\r^ -V/*)||| -h \\p-qh\\). 

By the triangle inequality, we obtain the estimate (14.281) . In addition, using the 
approximation property, we see that (14.291) holds. □ 

We can prove the error estimate of optimal order by the Aubin-Nitsche duality 
argument. 

Theorem 10 (L^-error estimate). Let the notation he the same as in Theorem\^ 
If {u,p) G Ff^+^(S1) X il^+^(n), then we have 

(4.34) ||ti -Uh\\< Ch^^'^ (|M|fj'=+2(n) + b|Hfc+i(n)) . 


Proof. We consider the following adjoint problem: find (il’.n) G (H^(Ll)r\H},(Q))x 
{H^in)nLl{n)) such that 

—Aif) + Vtt = u — Uh, 
divi/j = 0. 


Note that \'iP\h'^(q) + |7r|//i(Q) < C\\u — Uh\\. From the adjoint consistency (14.31) . 

the solution of the problem satisfies 

(4.35) 

ah{vh,Vh;i>,ip\rh) + bh{vh,Vh;Tr) = {u - Uh,Vh)n y{vh,Vh) £ x V,^, 

bh{ip,'ip\rh',(lh) =0 Mqh^Ql- 

Let (i/j/i, '0/1, TTh) G X X Qft be an approximation to (0,0|r;., tt) satisfying 

1(0 - 0/1,01 r;. -0/l)lll < C'/l|0|i^2(n), \\7T-TTh\\ < Ch\TT\Hl(n). 

Taking Vh = u — Uh and Vh = w|rh — Uh in (14.351) . 

\\u - UhW'^ = ah{u- Uh,u\r,, - ■u/i;0,0|rfe) + bh{u- Uh,u\T^ - Uh','^) 

= ah{u- Uh, u\r,^ -Uh;ip- if>h,ip\rh - ^h) 

Pbh{u- Uh, M|r^ -Uh',T^ - TT/i) 

< Cl\{u-Uh,u\r,, -M/i)|||(|||(0 -0,„0|r;. -0/i)lll + lk-7r,,||) 

< Ch^^'^ (|w|iT'^+2(n) + bli/fc+qn)) \\u - Uh\\. 

Thus we obtain the assertion. □ 
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5. Relations with the nonconforming Gauss-Legendre finite element 

METHOD 


5.1. The Gauss-Legendre element. The approximation space of the (fc + l)-th 
Gauss-Legendre element for velocity is defined by 

(5.1) = {Vh e = 0}, 

where [-J is a jump operator(see [2] for example). This space is known as the 
Crouzeix-Raviart (k = 0), Fortin-Soulie|19] (k = 1) or Crouzeix-Falk [T2](fc = 2) 
element. Note that is continuous at the (fc-l-l)-th order Gauss-Legendre 

points, and thereby PkVh is single-valued on F^. 

The nonconforming Gauss-Legendre finite element method reads as: find £ 

X such that 

(5.2a) (Vul,Vvh)rH-((iiv^h,p*h)rH = if,^h)n Vu/, e 

(5.2b) (divu);, qh)u =0 V% e Qt 

The nonconforming method is well-posed for k = 0,1,2. For fc > 3, it is the case 
under some assumption on a mesh, see [3j Lemma 3.1]. We here assume that the 
method (|5.2I1 is well-posed for simplicity. In our analysis, we will use the following 
inf-sup condition for the Gauss-Legendre element, see also [3]. 


Theorem 11 (Discrete inf-sup condition for the Gauss-Legendre element). There 
exists a constant f3 > 0 such that 


(5.3) 


sup 


{divVh,Ph)n 

\Vh\l,h 


> PkhW 


yqh e Ql 


5.2. A relation between the reduced-order HDG method with the lowest- 
order approximation and the Crouzeix-Raviart element. It is known that 
there are relations between HDG methods and the conforming or nonconforming 
finite element methods for Poisson’s equation. The hybrid part the solution of 
the embedded discontinuous Galerkin(EDG) methodpO] with the continuous Pi 
element is identical to the conforming finite element solution on F/i. In [23], it 
was proved that the hybrid part of the solution of the reduced method using the 
discontinuous Pi-Pq approximation coincides with the nonconforming Crouzeix- 
Raviart finite element method at the mid-points of edges. In this section, we 
discover a relation between the reduced method with the Pi-Pq/Pq approximation 
and the nonconforming Crouzeix-Raviart finite element method. 

Let Vj^ be the usual Crouzeix-Raviart finite element space. The Crouzeix-Raivart 
interpolation operator for scalar-valued functions, H]] : L^lTh) —>■ Vj^, is defined as 


(5.4) 


H’JjDds = 


vds 'ie £ Eh- 
,Vd)'^ £ L‘^{Th), we set 


For a vector-valued function v = (ui, 

Ulv := {Ulvi,--- ,U*hVdf 


We are now in a position to prove the relation. 

Theorem 12. Let {uh,Uh,Ph) £ Vfj;xV^ X Q® be the solution of the reduced method 
(13.111 with k = 0 and {u’^,pD £ Vh X Ql be the solution of the nonconforming 



ANALYSIS OF A REDUCED-ORDER HDG METHOD FOR THE STOKES EQUATIONS 11 


Crouzeix-Raviart finite element method. Then we have 

IllUh=ul, Ph=P*h- 

Proof. Since is a piecewise linear polynomial, it follows that, by the Green 

formula, 

= {VUlvj,,vniun)r,. 

Choosing = T^hVh in (13.fall and noting that vu — ^hTh = 0 at the mid-points 
of edges, we have 

(5.5) {VTllun,VTllvu)T, - {divUlvu,pn)n = S V^. 

The equation (I3.1bl) can be rewritten as 

-{(\:wUh,qh)Th - {^h - Uh,qhn)gTh = -{^h,qhn)gT^ 

(5.6) =-{Illuh,qhn)grh 

= -{divUluh,qhn)rh- 

From (15.51) and (15.61) . we obtain the following equations to determine Uh- 

(5.7) {VUlu^, Vnlv^)r, - {divUlv^,pn)r, = (/, Ulv^)n e 

(5.8) {d:wIVluh,qhn)TH =0 Vg?, G 

which is nothing but (|5.2I) in the case of fc = 1. Therefore we have = uf and 

Ph=P*h- □ 


Remark 13. From the definition of the Crouzeix-Raviart interpolation, we see that 

j Uhds = J U’^Uhds = j u’^ds, 

which implies Uh and u'^ are equal at the mid-point of e € Sh. 


5.3. The limit case as r tends to infinity. We will show that, as the stabiliza¬ 
tion parameter t tends to infinity, the solution of the reduced method converges to 
the nonconforming Gauss-Legendre finite element solution with rate 0(t“^). To 
do that, we first prove the following key lemma. 

Theorem 14. There exists a constant C > 0 such that, for any {vh,Vh) G x 

yk 

(5.9) inf l{vh-Wh,Vh-PkWh)lh<C\{vh,Vh)\y 

Proof. Let us define G(V)f^^) = {{vh,PkVh) ■ Vh G which is a closed 

subspace of x V^. We consider the quotient space x Vjf / 

where the standard quotient norm is given by 

II (ll/t) I’ll) II yfc+l X V''=/G(W+l) ~ lll('*’li ~'Wh,Vfi — Pfc'nj/j)|||/j. 

h hi \ h I Wh&vf+^ 

We shall prove that \{vh,Vk)\j is also a norm on x /G{V^^^). Suppose 

that |(i;^,n;i)|j = 0, then we readily have Vf, G and Vh = PkVh- Hence it 

follows that {vfi,Vh) G G{Vi^~^^), which implies that {vh,Vh) is the zero element 
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of X ' ). Since any two norms on a finite-dimensional space are 

equivalent to each other, we conclude the assertion. 

□ 


Let us denote the solution of the reduced method p.ll) with the 

stabilization parameter r. By using Lemma 1141 it can be proved that 
converges to zero as r —)• oo with rate 

Theorem 15. If t is sufficiently large, we have 


( 5 - 10 ) \{ul,ul)\^ < Ct 1 ||/||. 

Proof. Let Vh S be arbitrary. Since Pk{vh — PkVh) vanishes on L^,, we have 


l(“L«Dlj.T = \{ul-Vh,ul - 

< Wiul-Vhffil - PkVh)\lr. 


Therefore, 


(5.11) 


l(“L«Dlj.r < inf l\{ul-Wh,ul - PkWh)\lr. 


Taking {vhffih) = {vn, PkVh) in (|3T|) . we have 


(5.12) {yul,Vvh)r,, + {drvVh,ul - ul)on - (divt>/i,p;))rfc = {f,Vh)n. 

Subtracting this from (l.l.lall by taking {vh,Vh) = gives 

l(«L“Dlj> ={f^K-Vh)n - {yul,V{ul-Vh))TH 

(5.13) - {dnul,ul - uDoPh - {9niul - Vh),ul - uDon 

+ (div(M;; - Vh),pl)n + {ul - ul,pln)ar^. 

Here we note that the energy norm is stronger than the ffi norm(see [23] for the 
proof): 

(5.14) ||^;,|| < Cl\{zk,Zk)\U y{zn, zn) G x 

and that by the standard argument, we see that and are uniformly bounded 
in the energy and ffi norms, respectively, i.e. 


ul\\H + \\pl\\<C\\f\\ 
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where C is independent of r. The terms on the right-hand side of (15.1311 are bounded 
as follows: 

\{f.ul-Vh)a\ < \\f\\\\ul-Vh\\ 

< C\\f\\\liul-^h,ul - PkVh)\lh, 
\{yul,y{ul-Vh))Tj < l|V<||||V« 

< C'll/lllll« - - PkVh)\lh, 

\{dnUl,ul - uDorJ < C\\Vul\\\{ul,ul)\j 

< C'll/lll«>«Dlj: 

\{dr,iul-Vh),ul-ul)arJ < C\\V{ul - Vh)\\\iul,ul)\i 

< C\liul-Vh,ul - PkVh)\lh\{ul,ul)\j 
\{diviul-vh),pl)n\ < \\'^iK-Vh)\\\\pl\\ 

< c:i\{ul-Vh,ul - PkVh)\lh\\f\\ 

\{ul - ul,pln)9rj < C'l(«L“Dljlbhll 
<c|K,iiDljll/ll- 

From the above estimates, we obtain, by using Lemma [Ml 

l«.«Dlj-r < C'dI/ll + l«.«Dlj) • inf \liK-'Wh,ul - PkWh)l\h 
<t^*(ll/ll + IK,«Dlj)IK,«Dlj- 

Since we can assume r > C*, it finally follows that 

T — 

The proof is complete. □ 

In the following theorem, we prove that {ul^,pj^) converges to with the 

same rate as \{u'^,uj^)\j, namely 

Theorem 16. If t is sufficiently large, we have 

(5-15) \ul-ul\i^h + \\p*h-pl\\ <CT-^\\f\\. 

Proof Choosing Vh = Vh & Vh nnd = PkVh in (13.11) . we have 
(5.16a) 

(Vul, Vvh)n + {dnVh, ul - ul)dn - (divu?„pDr^ = if, Vh)n Vu?, e 

(5.16b) (div<, qhhn + iK “ <lhn)dTH =0 Vg?, e Qt 

Let us denote 'uf, — if, = and = p)) —Ph- Subtracting (15.161) 

from (ICT) leads to 

(5.17a) iV$,l,Vvh)Th - {dnVh,ul - uDar,, - {AwVh,6l)TH =0 Vu/, G 
(5.17b) (div|^, qh)TH + fK - K, Qhn)au =0 V% G Ql- 

Let Wh G be arbitrary, then we have 

(va, va)r. = (v« - wh), ^a)n + (V(i2/. - ui),ya)u 

=:I + II. 
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By (I5.17all with vt = u*j^ — Wh, we get 

I = (dn.(ul - Wh),ul - uDoTh + (div(M;; - Wh),Sl)u 

= {dnil.ul-uDar^ + {dr,{ul-Wh),ul-ul)aTH 

+ (div|^, 5l)r^ + (div« - wu), 5l)r^ 

='■ Ii + I 2 + I 3 . + li- 

These terms can be bounded as follows: 


\h\<c\eh\i,h\{ui,K)h, 

Ihl < C\ul - Wh\i,h\{ul,ul)\j, 

\h\ = \-{ul-ul,6ln)an\ (by (|5.17bD) 

<C\{ulul)\Ml 

\h\ < C\ul - Wh\i,h\K\\- 

Therefore, we have 


|/| < C{\a\i,h + \\Simul,ul)\^ + C(|(«L«Dlj + \\Sl\\)\ul - 

Taking the infimum with respect to Wh and using Lemma [TTl we have 


(5.18) 


\I\<C{\{ul,ul)\^ + \a\i,h + \\SU)\i<,iil)\y 


On the other hand, by the inf-sup condition for the Gauss-Legendre element, we 
have 


(5.19) 


f3\\Sl\\ < sup 
= sup 


{divvh,5l)r^ 

\Vh\l,h 

(y^l,\/vh)n - {dnVh,ul - ul)au 


I'l’hll.h 


< Kh\l,h + G|(m)), ItJOIj. 


Combining (I5.18|) with (j5.19l) gives us 

(5.20) |/| < c'(|«,M;;)|j -h \{ui,ui)\y 

The term //is estimated as follows by using Lemma (TT] 


(5.21) 


\II\ < inf \l{ul-Wh,ul - PkWh)\lh ■ \$.l\i,h 


<C\{ul,ul)\i\a\i,h- 

From (15.201) and 1)5.21|) . it follows that 

\ih\l,h < C” (|(m)), mDIj + l(“h> “DIj- 


By Young’s inequality, we get < CKit)),■u)))|j. Using (15.191) again, we have 

< C'l('*^L“Dlj' Theorem [T5l we obtain the assertion. □ 
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6. Numerical results 

As a test problem, we consider the case of = (0,1)^ and 

f{x,y) = (471^ sin(27ry), d-TT^ sin(27ra;)(—1 + 4 cos(27ry)))^. 

The source term is chosen so that the exact solution {u,p) is 

u(x,y) = (2 sin^(7ra;) sin(7ry), —2 sin(7ra:) sin^(7r?/))'^, 
p{x,y) = 4?! sin(27ra;) sin(27r2/). 

We carry out numerical computations to examine the convergence rates of the 
reduced method. We employ unstructured triangular meshes and the Pk+i-Pk/Pk 
approximations for fc = 0,1, 2. The convergence history is shown at Table[Tl and the 
convergence diagrams are displayed in Figure [1] From these results, we observe that 
the convergence orders are optimal in all the cases, which agrees with Theorems [9] 
andfTOl 


Table 1. Convergence history of the reduced method. 


Degree 

Mesh size 

\\U-Uh\\ 

||V(m- 

Uh)\\ 

\\P-Ph 

II 

k 

h 

Error 

Order 

Error 

Order 

Error 

Order 

0 

0.2634 

1.789E-01 

- 

3.672E+00 

- 

2.307E+00 

- 


0.1414 

4.938E-02 

2.07 

1.901E+00 

1.06 

1.249E+00 

0.99 


0.0701 

1.200E-02 

2.01 

9.340E-01 

1.01 

6.054E-01 

1.03 


0.0373 

2.907E-03 

2.25 

4.596E-01 

1.13 

2.925E-01 

1.15 

1 

0.2634 

1.038E-02 

- 

8.486E-01 

- 

4.520E-01 

- 


0.1414 

1.639E-03 

2.97 

2.290E-01 

2.11 

1.259E-01 

2.06 


0.0701 

1.815E-04 

3.13 

5.003E-02 

2.17 

2.768E-02 

2.16 


0.0373 

2.295E-05 

3.28 

1.254E-02 

2.19 

6.859E-03 

2.20 

2 

0.2634 

1.338E-03 

- 

1.439E-01 

- 

3.105E-02 

- 


0.1415 

8.656E-05 

4.40 

1.900E-02 

3.26 

3.920E-03 

3.33 


0.0701 

5.089E-06 

4.03 

2.298E-03 

3.01 

4.310E-04 

3.14 


0.0373 

3.168E-07 

4.40 

2.858E-04 

3.31 

4.956E-05 

3.43 
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